Floquet analysis of the modulated two-mode Bose-Hubbard model 
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We study the tunneling dynamics in a time-periodically modulated two-mode Bose-Hubbard 
model using Floquet theory. We consider situations where the system is in the self-trapping regime 
and either the tunneling amplitude, the interaction strength, or the energy difference between the 
modes is modulated. In the former two cases, the tunneling is enhanced in a wide range of modula- 
tion frequencies, while in the latter case the resonance is narrow. We explain this difference with the 
help of Floquet analysis. If the modulation amplitude is weak, the locations of the resonances can be 
found using the spectrum of the non-modulated Hamiltonian. Furthermore, we use Floquet analysis 
to explain the coherent destruction of tunneling (CDT) occurring in a large-amplitude modulated 
system. Finally, we present two ways to create a NOON state (a superposition of N particles in 
mode 1 with zero particles in mode 2 and vice versa). One is based on a coherent oscillation caused 
by detuning from a partial CDT. The other makes use of an adiabatic variation of the modula- 
tion frequency. This results in a Landau-Zener type of transition between the ground state and a 
NOON-like state. 

PACS numbers: 03.75. Lm, 33.80.Bo, 42.50. Dv, 74.50.+r 



I. INTRODUCTION 

Ultracold atomic gases are novel systems with a high 
degree of controllability, making them very useful in the 
studies on quantum phenomena. The possibility to con- 
trol the parameters during experiments is essential, for 
example, in quantum information processing (e.g., Refs. 
[J, 0) an d matter- wave interferometry (e.g., Refs. [j- 
|5[). In this paper, we discuss the dynamics of an ultra- 
cold bosonic gas trapped in a time-periodically modu- 
lated potential. The dynamics of periodically modulated 
quantu m sy stems has attracted both theoretical (e.g., 
Refs. @-|32) and experimental (e.g., Refs. |33l - [4l| ) in- 
terest during recent years. It is known that a modulated 
system has resonances at which the tunneling is cither 
suppressed or enhanced. In the neighborhood of a reso- 
nance, the behavior of the system is very sensitive to the 
modulation frequency. 

The suppression of tunneling by modulating the energy 
difference between the modes is known as the coherent 
destruction of tunneling (CDT) 0-Q. CDT was discov- 
ered in Ref. [6|, where the motion of a charged particle 
in a lattice under the influence of an oscillating electric 
field was studied. It was shown that an initially localized 
particle remains localized in a one-dimensional lattice if 
the amplitude and frequency of the electric field are cho- 
sen suitably. In Ref. 0, CDT was found to occur in 
systems consisting of a particle subjected to a periodic 
force and trapped in a double-well potential. Recently, 
the coherent destruction of tunneling has been actively 
studied in the context of ultracold bosonic atoms (e.g., 

Refs. raEHIilliiliiHSIsiTdjI). 



Unlike the CDT, which is typically observed under the 



condition that the tunneling coupling is larger than or 
comparable to the interaction energy, the enhancement 
of tunneling by modulation can take place in a system 
where the interaction energy dominates over the tunnel- 
ing coupling. In the absence of modulation, the large in- 
teraction energy suppresses tunneling for energetic rea- 
sons [42, |43(. This leads to a very long tunneling pe- 
riod (the time needed for N particles to tunnel from one 
mode to another and back). However, by modulating 
the tunneling matrix element, it is possible to enhance 
the many-particle tunneling and thereby reduce the tun- 
neling period [28J]. In this paper, we analyze the reasons 
behind the enhanced tunneling with the help of a de- 
tailed Floquet analysis. In order to make the analysis 
more complete, we consider also systems where either 
the interaction strength or the energy difference between 
the modes is modulated. We find that these two meth- 
ods provide an alternative way to enhance tunneling. It 
is shown that the width of the resonance, that is, the 
range of modulation frequencies corresponding to the en- 
hanced tunneling, depends strongly on whether the tun- 
neling matrix element, the interaction strength, or the 
energy difference between the two modes is modulated; 
the resonance is wider in the former two cases. We ex- 
plain this difference with the help of Floquet theory and 
the eigenvalues of the non-modulated Hamiltonian. 

We analyze also the coherent destruction of tunnel- 
ing using Floquet theory It is known that CDT can 
be caused by modulating either the energy difference be- 
tween the modes or the interaction strength. We study 
only the Floquet spectrum of the former system because 
it has not received much attention in the literature, unlike 
the Floquet spectrum of the interaction-modulated sys- 
tem [22|, |26| . In addition to this, we present two ways to 



generate NOON-likc (Schrodingcr's-cat-like) states. The 
first is based on the CDT induced by a large-amplitude 
modulation of interaction strength, whereas the second 
makes use of a small-amplitude modulation of the tun- 
neling coupling. 

This paper is organized as follows. In SecHU we define 
the modulated Bose Hubbard Hamiltonian. In Sec. IIII1 
we give a short summary of the Floquet theory used in 
this article. Section [TV] discusses in depth the results 
of the Floquet analysis for systems in the self-trapping 
regime subjected to a small-amplitude modulation. In 
Scc.|Vj the coherent destruction of tunneling is examined 
using Floquet theory. It is also shown that NOON states 
can be created with the help of partial CDT. In Sec. lVIl a 
way to create NOON states using adiabatic sweep across 
an avoided crossing is presented. Finally, the conclusions 
are in Sec. IVIII 



and interaction strength are defined as 

U(t) = U + U 1 sm(wt + <f> u ), 
V(t) =V + Vism(wt + <f> v ), 



(7) 
(8) 



where Uq,Vq (Ui,Vi) are the amplitudes of the static 
(time-dependent) part of the interaction and the tilt, re- 
spectively. In the above equations, uj is the modulation 
frequency and 4>j,4>u, and 4>v ar e the phase offsets. In 
this paper, time is measured in units of 



„ 71" 

Jo 



(9) 



which is the tunneling period in the absence of the inter- 
action (U = U\ = 0) and the tilt (V = Vi = 0). Here 
and in what follows, we set H = 1. 



II. TIME-PERIODICALLY MODULATED 
TWO-SITE BOSE HUBBARD HAMILTONIAN 

We consider a system described by the two-mode Bose- 
Hubbard Hamiltonian. For definiteness, we assume that 
this model is realized physically by bosons trapped in 
a double-well potential. We consider situations where 
cither the tunneling amplitude, the interaction strength, 
or the energy difference between the wells is modulated 
periodically in time. This system is described by the 
Hamiltonian 



H(t) = -J(i)(c t 1 c 2 + 4ci) + ^(c t 1 c t 1 fiici + 

+ -^-(4ci-clc 2 ) 

= -2J(t)S x + U{t)S 2 z +V(t)S z . 



C2C2C2C2) 

(1) 

(2) 



Here J is the tunneling matrix element, U is the on-site 
interaction, and V is the energy difference between the 
wells (tilt). We have introduced the SU(2) generators 
defined as 



b x - 


= ^i^+^ci) » 


(3) 


Sy = 


= Yi^ 1&2 ~ £ 2Cl) , 


(4) 


s, -- 


= o(clci -c\c 2 ) , 


(5) 



where Ci{c[) annihilates (creates) an atom in mode i. 

We define the time-dependent tunneling matrix ele- 
ment as 



J(t) = J [1 + A j sin (ut + <f> j)], 



(G) 



where Jo is the amplitude of the time-independent part 
and Aj <S [0, 1] gives the relative amplitude of the time- 
dependent tunneling matrix clement. The modulated tilt 



III. FLOQUET OPERATOR 

If the Hamiltonian H is periodic in time, Floquet the- 
ory provides a powerful tool to analyze the dynamics of 
the system. In the following, we denote the modulation 
period of H by T u . In our case, the modulation is sinu- 
soidal and hence T u = 2tt/oj. According to the Floquet 
theorem (see, e.g., Ref. [4J]), the time-evolution opera- 
tor Uh determined by the Hamiltonian of Eq. $Z§ can be 
written as 



U k (t) = M(t) 



-UK 



(10) 




where M is a periodic matrix with minimum period T u 
and M(0) = I and K is a time-independent operator. We 
define the Floquet operator F as 

(11) 
1 l H(t)dt \ , (12) 

where T is the time-ordering operator. At times t = nT u , 
where n is an integer, we get [/^(nT^) = e -" lT ^ K — p n . 
The Floquet operator is a mapping between the state at 
t = and the state after one modulation period T u = 
2tt/w: <&(T u ) = F^(0). The columns of the Floquet op- 
erator F can be obtained by following the time evolution 
of the basis states for one modulation period. Each time- 
evolved basis state forms a column of the matrix F. The 
Hilbcrt space of a two-mode system containing N bosons 
is C N+1 . The basis of this Hilbert space can be chosen to 
be {| AN) ; A7V = -N, -N + 2, -N + 4, . . . , N}, where 
|AJV) is a state with (N + AN)/ 2 particles in mode 1 
and (N — AN)/2 particles in mode 2. Any pure state of 
the system can be written as 



N 



■>!• 



J2 C AN \AN), 



(13) 



AN=~N 



where the amplitudes {Can} are complex numbers. If 
N is even (odd), AN takes only even (odd) values. 

In order to characterize the eigenstates of the Floquct 
operator, we define the parity operator P as 



P\AN) 



AN). 



(14) 



It can alternatively be written as P = (— i) N e mSx . The 
eigenvalues of P are 1 and — 1, corresponding to even 
and odd parity, respectively Because P^S Z P = —S z , 
the Hamiltonian, and consequently the time-evolution 
operator, commutes with P if the tilt vanishes. Then 
the eigenstates of F are also eigenstates of P and either 
Can = C-an or Can = -C_aw holds for the compo- 
nents of the eigenvectors of F. In the former case, the 
cigenstate has even parity and is said to be symmetric, 
while in the latter case the parity is odd and the eigen- 
state is called antisymmetric. Furthermore, the absolute 
values of the coefficients {Can} of an eigenstate have 
maxima at AiV = ±fc, where k > is an integer. We 
denote such an eigenstate by ip k , where + (— ) indi- 
cates that the eigenvector is symmetric (antisymmetric) . 
If Jo <C UqN and the signs of the eigenvectors are defined 
appropriately, we get 



*$> 



T2 m± 



N)), 



(15) 



which is valid to zeroth order in Jq/I/qN. For non-zero 
Vq or Vi , the Floquet eigenstates are neither exactly sym- 
metric nor antisymmetric because S z is not invariant un- 
der the parity operator. However, since the time average 
of the S z term is zero (we assume that Vq = 0), the Flo- 
quet eigenstates are almost symmetric or antisymmetric, 
provided V\ is small. We thus use the notation ip\ also 
in this case. Note that, in the case of large-amplitude 
modulation of the tilt, the Floquct eigenstates cannot be 
classified in this way. 

The Floquet operator is a unitary operator, and there- 
fore the eigenvalue corresponding to the eigenvector ipS 

. ,(±) 
can be written as e l>Pi . The eigenvalue equation be- 
comes 



In this paper, we call 
quet eigenvalue. 



(±) 



o *0. 



(±) 







(±) 



(16) 



a(±) 



G [— vr, 7r) the phase of a Flo- 



Assume that the initial state is \N) ra (-0 
V>jy )/-\/2. At t = riT^ (n G N), the state reads 



(+) 

N 



F n \N) 



in<j>\ 



v/2 



(■</■ 



( + ) , An[<fc 



N 



&w 



(17) 



ifH<4 -) 



6^| w it, we get F n \N) ra | -TV); that is, the 
system has tunneled from \N) to | — N). In this paper, 
we define the tunneling period as the time needed for the 
system to tunnel from \N) to | — N) and back. In terms 



of the phases of the Floquet eigenvalues, the tunneling 
period reads 



T 



2tt7L 



4" } 



4 +) l 



(18) 



Increasing \4> N — 

vice versa. When 
verges. 



i> N | reduces the tunneling period and 
(f) N = <f) N , the tunneling period di- 



IV. TUNNELING PERIOD AND FLOQUET 
ANALYSIS IN THE SELF-TRAPPING REGIME 

In this section, we consider the tunneling of bosons in 
the self-trapping regime characterized by UoN/2J ^S> 1. 
Assume that in the initial state all N particles are either 
in site 1 or site 2. The reduction of the interaction en- 
ergy by single-particle tunneling is of order ~ UqN . This 
reduction cannot be compensated by the increase of the 
kinetic energy, which is approximately given by ~ Jo- 
As a consequence, single-particle tunneling is suppressed 
(self-trapping), and all N particles stay in the same well 
for a long time. In this situation, oscillations between 
the states IN) and | — N) occur through higher-order co- 
tunneling [45( . In Ref. [28[ it was found that the tunnel- 
ing period of the higher-order co-tunneling can be dras- 
tically reduced by modulating the tunneling matrix ele- 
ment J. As we show here, a similar phenomenon can be 
seen when the tilt is modulated (we set Vq = 0). In Figs. 
Q] and [H we show the tunneling period T for the modu- 
lated tunneling matrix element and tilt, respectively. As 
an example, we have set N = 5 and Uq/Jq = 4 in the 
both cases. We see that the behavior of J as a function of 
the modulation frequency u> depends strongly on whether 
J or V is modulated. This difference can be explained 
using Floquct analysis. 

Before analyzing the system in detail, we first summa- 
rize two key points. One is the parity of the operator 
whose coefficient is modulated, and the other is the shift 
in the phases of the Floquet eigenvalues due to an avoided 
crossing. The parity of S x is even and that of S z is odd. 
Therefore, S x couples Floquct eigenstates of the same 
parity, and S z couples those of the opposite parity. This 
means that in the case of modulated J [V], there is an 
avoided crossing between Floquet eigenstates of the same 
[opposite] parity. The differences in the behavior of the 
tunneling period can be attributed to the parities of the 
states undergoing an avoided crossing. Below we show 
that usually <j^ ] > 4~ ] (^~ ] > 4 +) ) holds for odd 
(even) N. However, this is not necessarily the case near 

avoided crossings where the values of (f>\ are shifted. 
These shifts lead to either suppression or enhancement 
of the tunneling. 

We have chosen TV = 5 and Uo/Jo = 4 in this section. 
The results can, however, be straightforwardly general- 
ized to any value of N and Jq/Uq, as long as N > 1 
and U N/J Q > 1. Wc remark that Figs. Q] 13(b), and 
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FIG. 1: (Color online) Tunneling period T in the case of mod- 
ulated tunneling matrix element J for N — 5, Uo/Jo — 4, and 
Aj = 0.1 (and V = Vi = U\ = 0). We have set <jjj = 0, 
but there are no noticeable differences for different values of 
4>.j. There is a drastic reduction of T in a wide range around 
uj/ Jo ~ 16. Very narrow resonances in the region cj/Jq < 10 
are not shown. The vertical red dotted lines and arrows show 
the positions of the resonances evaluated from Eq. ()22|) using 
the energy eigenvalues of the time-independent Hamiltonian. 
This figure is adopted from Ref. [281 ] - 



A. Time-independent Hamiltonian 

If the modulation amplitude is small and the system is 
not near an avoided crossing, the Floquet eigenstates and 
eigenvalues turn out to be close to the ones determined 
by the time-independent part of the Hamiltonian, given 
by 



Ho 



-2JqS x + UqS z . 



(19) 



As a consequence, some important properties of the mod- 
ulated system, such as the positions of the resonances, 
can be explained by analyzing the spectrum of Hq . 

We assume that UqN S> Jo and Vq = 0. In order 
to compare the time-evolution operator of the original 
time-dependent modulated system with that determined 
by Hamiltonian (fl~9| . we define the Floquet operator Fo 
corresponding to Hq as 



F n 



1 J- U) — 



(20) 



We denote the eigenvalues of the time-independent 
Hamiltonian by Ej\. k , where we use the same indexing 
as in the case of the eigenvectors of the Floquet operator 
F. The phases of the Floquet eigenvalues are given by 



b (±) 

y 0;fc 



P (±) T 



mod 2-7T. 



(21) 




FIG. 2: (Color online) Tunneling period T in the case of mod- 
ulated tilt for N = 5, Uo/Jo = 4, and Vi/Jo - 0.2 (and 
Vo = A j — Ui =0). We have set <f>v — 0, but there are no 
noticeable differences for different values of cj>y. The vertical 
red dotted lines and arrows show the positions of the reso- 
nances evaluated from Eq. (|22[1 using the energy eigenvalues 
of the time-independent Hamiltonian. 



the top panel of Fig. 2] are taken from Ref. [28( , but the 
detailed Floquet analysis of the J modulation, as well as 
the entire analysis of the effects of the tilt and interaction 
modulation, has not been presented elsewhere. 



We find that E$ < E$ for odd N and E$ > E$ 
for even N. Because of the minus sign in Eq. (|2"Tj) , the 
opposite holds for the phases of the Floquet eigenvalues 
(f> . k . The situation is similar in the presence of a small- 
amplitude modulation, and thus, normally, <f>\ > (fy. 



«.(-) 



> 



d+h 



for odd (even) N. Now E . k > E 0; i if fc > 
I > 0. Using this and the equation 9 w </> . fc = E Q . k 2tt/uj 2 , 

we see that </> ofc , anc ^ therefore also <p k , increases faster 
as a function of to the larger k is. This means that if k > I, 
4>k approaches <f>i from below as to increases [see Fig. 131b) 

for an example]. The phases {<f> k } cross repeatedly as 
ui increases. A crossing occurs when ui satisfies 

nwwI^ojfc-^Djjl, (22) 

with n — 1,2,3,.... The last crossing between <pk and 
(f>i is at u ~ \Ek — Ei\. In the limit of uj — > oo, the 
phases of all the Floquet eigenvalues approach zero from 
the negative side. 

In the specific case N = 5 and Uo/Jq = 4, corre- 
sponding to Figs. Q] and O the crossing points between <f>^ 
and the other phases are at uj/Jq = 28.01 (crossing with 



^[ +) ), 22.63 



b (-h 



15.93 
are obtained from Eq. 



d+h 



and 15.31 



with n = 1. Note that E, 



4 } )- These 

(+) 



and E§ , and thus 



/.(+) 



and 



are almost identical. 



The crossing points corresponding to n = 1 and 2 in 
the region u>/ Jo > 10 are shown by the vertical red dot- 
ted lines and arrows in Figs. Q] and f2J We see that the 
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FIG. 3: (Color online) Results of the Floquet analysis for the case of Fig. [T] i.e., modulated hopping parameter J. Panel (a), 
which is the same as Fig. Q] (apart from the range of the horizontal axis), shows the tunneling period T as a function of the 
modulation frequency uj. Panel (b) shows the phases of the Floquet eigenvalues, 0|aat|, as a function of ui. Panel (c) shows 
the schematic behavior of the Floquet eigenvalues near the crossing points in panel (b). Each resonance observed in panel (a) 
corresponds to one of the three types of crossings shown in panel (c). In panels (a) and (b), crossings of types 1, 2, and 3 are 
labeled by the magenta solid, green dashed, and blue dotted curves, respectively. 



positions of all the resonances shown in Figs. [TJ and [H 
are well explained by the energy eigenvalues of the time- 
independent Hamiltonian. Based on this fact, we can say 
that the positions of the crossings are the same irrespec- 
tive of the modulated variable. 



tcm, there are N — 2 different types of crossings between 

s j46(. Because now N = 5, we 



j(±) 



■N 



and other 



d±) 



have three types of crossings; each resonance corresponds 
to one of these. In the following, we analyze in detail each 
of these three crossing types. 



B. Modulated J 

Next we consider the modulation of the tunneling ma- 
trix element; see Figs. Q] and [3] The value (ipi\S x \ipj) can 
be non-zero only if the Floquet eigenstates ipi and tpj 
have the same parity. Consequently, there is an avoided 
crossing between eigenstates with the same parity. 

In Fig. [3jb), we show the phases of the Floquet eigen- 
values as a function of ui/ Jq for the parameters used in 
Fig. Q] [and Fig. EJa)]. From Fig. [3] we see that a large 

change of T occurs when cp^ ' crosses the other <t>\ AN < 's 
[circles and ellipses in Figs.EJa) and^b)]. The behavior 
of the phases of the Floquet eigenvalues near the crossings 
is schematically shown in Fig. G^c). In an iV-particle sys- 



1. Type 1 

A crossing between (jr 5 and </>g leads to a reduction 
of T in a wide range of u> around lo/Jq — 16. This crossing 
is indicated in Figs. G2a) andEJb) by the solid magenta 
circle. The detailed structure of the crossing is shown 
schematically in the top figure in Fig. j3jc) . Since (j> 3 
and </>3 are almost equal, the avoided crossings between 
(f> 3 and <f> 5 and between cp^ and </>g occur almost 
simultaneously [the red solid circles in Fig. HJc)]. Be- 
cause 4>\ > 4>\ for odd N outside the crossing region 
(see Sec. IIV Ap and S x couples Floquet eigenstates with 
the same parity, the first avoided crossing occurs between 



(j) 3 ' and (j) 5 [the left red solid circle] as the modulation 
frequency increases. Due to the repulsion between these 
two levels, the splitting between <p§ ' is increased near 
the avoided crossing and thus the tunneling period is re- 
duced. The second avoided crossing takes place between 
(f>\ ' and 4> 5 ' [the right red solid circle]. Note that, af- 
ter the first avoided crossing, the states ip% and ^5 
have been interchanged [between the red solid circles] and 
the energy splitting between cf> 5 remains large until the 
second avoided crossing at which ^3 and %p§ ' are in- 
terchanged. These successive avoided crossings lead to 
a reduction of the tunneling period in a wide range of 
uj/Jq. 



proach the crossing from the small uj/Jq side. In this 
case, we have an avoided crossing between <j)\ and <fi 5 . 

The phase ^5 , which is located above <p§ far from the 
crossing, is pushed downward due to the avoided cross- 
ing with (/>[ , and thus <j) 5 and cf> 5 cross each other 
(the small blue solid circle). This leads to the divergence 
of the tunneling period. After this, there is an avoided 
crossing between 4>\ ' and (jr 5 (the large red solid cir- 
cle), leading to an enhancement of the quasienergy split- 
ting between 4> b . This yields a reduction in the tunnel- 
ing period. After the avoided crossing, the states ip[ 
and ?/>g are interchanged. 



2. Type 2 



4. Type 1': Small A j 



Because of the large quasienergy splitting between <j>\ 
and 4>i , the points where 4> b crosses 4>{ and <f>\ ' are 
far apart. We call a crossing between <j) 5 and <£>[ a, 

type 2 crossing and that between </>g ' and <j>\ a type 3 
crossing. With increasing w, a type 2 crossing first yields 
a reduction and then an enhancement of the tunneling 
period. We show the schematic structure of a type 2 
crossing in the middle figure in Fig.[3Jc). The resonances 
around uj/Jq ~ 11 and uj/Jq ~ 23, indicated by the green 
dashed curves in Figs.[3ja) and^b), correspond to type 
2 crossings. 

Suppose that the crossing is approached from the small 



k(+) 



> 



fc(-) 



as 



uj/Jq side. Far from the avoided crossing 
explained in Sec. IIVAI Since S x couples Floquet eigen- 
states with the same parity, <j>\ undergoes an avoided 
crossing with <p§ (the large red solid circle). Near the 

avoided crossing, the energy splitting between <f> 5 is en- 
hanced, which leads to the reduction of the tunneling pe- 
riod. Just after the avoided crossing (around the vertical 
dashed line), the states ip[ and ip 5 are interchanged 
and, unlike in the usual situation, <j> 



(-) 



> 



k(+) 



Since 



4> 5 is larger than cf> 5 ' far from the avoided crossing 

also on the large uj/Jq side, <f> 5 and (j> 5 cross each 
other (the small blue solid circle), which yields a diver- 
gence of the tunneling period. 



In the top panel of Fig. QJ we show the tunneling pe- 
riod T for various values of the modulation amplitude Aj. 
With decreasing Aj, the resonance around uj/Jq ~ 16 be- 
comes narrower and finally separates into two resonances 
(see the case Aj = 0.01 shown by the blue dashed-dottcd 
line). The schematic behavior of the phases of the Flo- 
quet eigenvalues near the crossing is shown in the bottom 
panel of Fig.HJ We call this a type 1' crossing. The major 
difference between type 1' and type 1 crossings is the ex- 
istence of two points where cj) 5 cross each other. These 
are indicated by the small blue solid circles, and they are 
located between the two avoided crossings (the large red 
solid circles). One can also view a type 1' crossing as a 
combination of type 2 and type 3 crossings. 

When Aj is small, the coupling between the two states 
that undergo an avoided crossing is small. Thus the dif- 



ference 



a(+) 



remains very small even near the 



avoided crossing. Therefore, unlike in a type 1 cross- 
ing, the inverted configuration of <f> 5 (i.e., the situation 

4> 5 > </>5 ) cannot be sustained throughout the region 
between the two avoided crossings. This leads to the ap- 
pearance of two crossing points, indicating a diverging 
tunneling period. In Ref. [281 ] it has been shown that the 
divergences are present if Aj < N^ 1 (Jq/Uq) n ^ 3 (N — 
l)(AT-2)/(JV-3)!. 



Modulated V 



3. Type 3 

As opposed to a type 2 crossing, a type 3 crossing 
(crossings between (j> 5 and <j>\ ) gives first an enhance- 
ment and then a reduction of the tunneling period. The 
resonances at uj/Jq ~ 14 and uj/Jq ~ 28, which are indi- 
cated by the blue dotted ellipses and circles in Figs. [3Ja) 
and EJb), correspond to type 3 crossings. A detailed 
schematic structure of a type 3 crossing is shown in the 
bottom figure in Fig. EJc) . Suppose again that we ap- 



As can be seen from Figs. \T\ and [2 the tunneling pe- 
riod behaves differently when V is modulated. In Fig. [51 
a noticeable change in the tunneling period T can be 
seen around uj/Jq ~ 16. There are also small, narrow 
resonances at uj/Jq ~ 22.5 and uj/Jq ~ 28 [47|. Unlike 
in the case of modulated J, the resonance at uj/Jq ~ 16 
is not a wide and smooth reduction of T for any value of 
the modulation amplitude V\. 

As in the case of modulated J, the resonance around 
uj/Jq ~ 16 is caused by a crossing between <f>\ and 
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FIG. 4: (Color online) The top panel (taken from Ref. |2g]) 
shows the tunneling period T in the case of modulated tun- 
neling matrix element J for various values of the modulation 
amplitude. The amplitudes used are Aj = 0.5 (red dotted 
line), 0.1 (black solid line), 0.05 (green dashed line), and 0.01 
(blue dashed-dotted line). The other parameters are the same 
as in Figs. [QandE] TV = 5 and U /J = 4 (and Vo = Vi = 0). 
The tunneling period does not depend on <f>j noticeably, and 
here we have set cj>j = for definiteness. The bottom panel 
shows the schematic behavior of the phases of the Floquet 



a<±) 



/,(±) 



eigenvalues near the crossing between <p 5 and (fi^ for small 
values of Aj (corresponding to, e.g., Aj = 0.01 in the top 
panel) compared to the type 1 case shown in Fig. |3jc) . We 
call this a type 1' crossing. 



</>3 \ However, unlike S x , the operator S z has odd par- 
ity, and it thus couples Floquet eigenstates of opposite 
parity. In Fig. [SJ we show the schematic behavior of the 
Floquet eigenstates near uj/Jq ~ 16. Suppose that the 
crossing is approached from the small uj/Jq side. As uj/Jq 
increases, the states ips an d V's undergo an avoided 



crossing, and 



d+) 



is pushed downward. Far from the 



d+) 



> 



/,(-) 



holds. Because 



avoided crossing the relation 

of this and the fact that a^ ' is pushed downward, the 

phases 4> 5 and </>5 cross (the left small blue circle) be- 
fore the avoided crossing (the left large red circle). After 
the first avoided crossing, the states V4 an d %p§ are 



*v 




FIG. 5: (Color online) Schematic behavior of the phases of 
the Floquet eigenvalues near the resonance around u)/ Jo — 16 
in the case of modulated V shown in Fig. [2] 



crossing (the right large red circle), and the correspond- 
ing states arc interchanged. Because now cpf' > <fy 5 , 
these phases cross after the second avoided crossing (the 



/.(+) 



> 



4~> far away 



right small blue circle), so that 
from the crossing. The two crossing points and the two 
avoided crossings correspond to the two divergences and 
the two reductions of the tunneling period, respectively. 
These are shown in Fig.[2]near uj/Jq ~ 16. Because in the 
present case the avoided crossings occur between Floquet 
eigenstates of opposite parity, the phases </4 ' necessar- 
ily cross each other outside the region of the successive 
avoided crossings. For this reason, a smooth reduction of 
T in a wide region of the modulation frequency ui cannot 
be achieved by modulating the tilt. This is one of the 
major findings of this paper. 

We note that all the other resonances are also much 
narrower than in the case of modulated J. This is be- 
cause the operator S z , which is related to the tilt, docs 
not contribute to the single-particle tunneling, unlike 
S x . The range of uj characterizing the width of the res- 
onance is comparable to the quasienergy separation at 
the avoided crossing. This is approximately proportional 
to | (■05 I^IV'IS)! 2 m the case of J modulation and to 
| ("05 |'S z |V'L4g)| 2 in the case of V modulation. The latter 
is smaller than the former by a factor ~ (Jq/Uq) 2 . This 
will be discussed in more detail in Sec. IVII 



interchanged. Next, 



d+) 



and 



/,(-) 



undergo an avoided 



D. Modulated U 

Finally, we consider the case in which the on-site in- 
teraction strength U is modulated weakly (JJ\/Uq <C 1). 
The Hamiltonian in this case is H = —2JqS x + U(t)S 2 , 
with U(t) given by Eq. ([7]). Since this Hamiltonian can 
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FIG. 6: (Color online) Tunneling period T in the case of U 
modulation (red solid line). For comparison, the tunneling 
period corresponding to J modulation is also shown (blue 
dashed line). Here N — 5, Uo/Jo = 4, and Vo = V\ = 0. In 
the case of J modulation j4j = 0.1 and Ui = 0, while in the 
case of 17 modulation Aj = and Ui/Jo = 0.4. 



FIG. 7: (Color online) Tunneling period in the weak interac- 
tion regime with large-amplitude modulation of the tilt. The 
parameters are N = 5, Uo/Jo = 0.1, Vo = 0, and Vi/Jo = 10 
(and A j — Ui — 0). We have set 0y = in this calcula- 
tion, but the behavior of T does not depend noticeably on 
<j>v- The vertical red dotted lines correspond to the values of 
lo /Jo which give Jo{Vi/uj) — 0. 



be rewritten as 

H{t) = A(t) [-2J eS (t)S x + U S 2 Z 
with A(t) = 1 + (Ui/Ua) sin (ut + 4>u) and 



Je&(t) — Jo 



1 + — sin (u)t + 4>u + ti") 
Uo 



(23) 



(24) 



we can expect that the dynamics can be reproduced by 
modulating J with the amplitude Aj = U±/Uo and phase 
4>J = 4>u + f instead of modulating £/. 

This observation is confirmed by the result shown in 
Fig. [51 where we compare the tunneling period T as a 
function of u) in the cases of modulated J and modulated 
U. In this example N = 5 and Uo/Jo = 4. The result 
for the modulated J is taken from Fig. UKAj = 0.1). By 
setting U\ = AjU , i.e., C/i/J = AjU /Jo = 0.4 in the 
present case, these two results almost coincide with each 
other. Note that since T does not noticeably depend on 
the phase of the modulation, U\ = AjUo is a sufficient 
condition for the behaviors of the tunneling periods to 
coincide. 

An analysis of the phases of the Floquet eigenvalues in 
the case of the U modulation shows that the schematic 
behavior of these phases around each resonance is the 
same as in the case of the J modulation shown in Fig. 
OJc). This can be understood by noting that S x and S 2 Z 
have the same parity. 



V. COHERENT DESTRUCTION OF 
TUNNELING AND THE FLOQUET SPECTRUM 

In this section, we first study a system characterized 
by a weak interaction and a large-amplitude tilt mod- 
ulation, concentrating on the properties of the Floquet 



spectrum. After this we examine the effects of a large- 
amplitude modulation of the interaction. The Floquet 
spectrum of this system has been analyzed elsewhere (see 
Refs. |22|,[26[) and will not be discussed here. Instead, we 
propose a way to create NOON states using selective tun- 
neling originating from the modulation of the interaction 
strength. 



A. Modulated V 

Next, wc consider a case where the interaction is weak, 
UN/ Jq < 1, and the amplitude of the modulation of the 
tilt is large, V\/Jq ^> 1. Wc assume that the tunnel- 
ing matrix element J and the interaction strength U are 
time- independent, that is, Aj = and U\ = 0. Further- 
more, we set Vo = 0. In this case, it is well-known that 
the effect of the modulation of the tilt can be approxi- 
mately described by a rcnormalized tunneling term. In 
more detail, the original tunneling term T = —2JS X is 
replaced by an effective one @, 0, dl, 03, HJ 0, HI : 



T, 



eff 



2Jo 



Mm 



— cos 4>v 



b.r 



COS <pv 

u 



s, 



■■}• 



(25) 



where Jo is the zeroth-order Bessel function (see Ap- 
pendix [A] for the derivation). Coherent destruction of 
tunneling takes place when Vi/ui is equal to one of the 
zeros of Jo- In the rest of this section, we discuss CDT in 
terms of the Floquet eigenvalues. This discussion holds 
for any N > 1. 
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co// 







Floquct cigcnstatcs that do not have maximum ampli- 
tude at the AN = component but that always have 
a non-zero AN = component, and 3) N/2 Floquct 
eigenstates that do not have maximum amplitude at the 
AN = component and where this component becomes 
zero when CDT occurs. The trio consists of all the three 
types, and the (N — 2)/2 pairs consist of the second and 
third types. We note that, for even N, the degeneracies 
of the trio and of all the pairs are incomplete provided 
Uq 7^ [48[, while all the pairwisc degeneracies are com- 
plete for odd N. Consequently, CDT is more complete 
for odd N than even N. 



Modulated U 



FIG. 8: (Color online) Phases of the Floquet eigenvalues in 
the case of Fig. [Jj The parameters are N — 5, Uo/Jo = 
0.1, V = 0, and Vi/J = 10 (and Aj = Ui = 0). In this 
calculation we have set 4>v = 0, but the phases of the Floquet 
eigenvalues do not depend on 4>v- The vertical red dotted 
lines correspond to the values of to /Jo which give Jo(Vi/u)) — 
0. The red arrows show the actual positions of the peaks of 
T (see Fig. 0). 



In Fig. we show the tunneling period T as a func- 
tion of the modulation frequency uj in the regime of weak 
interaction and large-amplitude modulation. In this cal- 
culation, we have set TV = 5, Uq/Jo = 0.1, Vq = 0, 
Vi/Jo = 10, and aV = 0, and in the initial state all 
particles are in site 1. The first five zeros of Jo(Vi/ui) 
are at Vi/u = 2.40, 5.52, 8.65, 11.79, and 14.93: they 
correspond to uj/J = 4.16, 1.81, 1.16, 0.848, and 0.670, 
respectively. These frequencies are shown by the vertical 
red dotted lines in Fig. [7J There is good agreement be- 
tween these dotted lines and the actual positions of the 
peaks of T. 

In Fig. [5J we plot the phases of the Floquet eigenvalues 
for the parameters used in Fig. [7] When the CDT occurs, 
the phases gather in pairs, the phases in each pair being 
almost equal, and all the pairs gather in a narrow region 
(red arrows in Fig. [5]). This behavior can be understood 
by noting that the Hamiltonian is effectively ~ UqS% at 
the point where CDT occurs, and thus A7V becomes a 
good quantum number, with a twofold degeneracy with 
respect to ±AN. 

Finally, we discuss the difference between even and odd 
N cases. For even N, the number of the Floquet eigen- 
values is N + 1, which is odd. Therefore, when CDT 
occurs, the Floquet eigenvalues are grouped into one trio 
and (N - 2)/2 pairs [cf. (N + l)/2 pairs for odd N}. 
A key point is that, for even N, there is a Fock state 
\AN = 0), which does not have a degenerate pair, unlike 
the other Fock states. In this case, the Floquet eigen- 
states near the value of ui at which CDT occurs can be 
classified into three types: 1) one Floquet eigenstate that 
has maximum amplitude at AiV = component, 2) N/2 



Due to the non-linear dependence of the interaction on 
AN, the CDT caused by a large-amplitude modulation of 
the interaction strength (Ui ^> Jo, Uq) is state dependent 
[2a ]. Here we assume Aj = V = for simplicity. In 
this case, a condition for partial CDT between the states 
\AN = to) and \AN = m — 2) (to. is a positive integer) 
is 



Jo 



— (m-1) 



= 0; 



(26) 



see Appendix iBl for the derivation. 

Unlike in the case of modulated V shown in Fig. [51 only 
the Floquct eigenstates relevant to partial CDT show the 
degeneracy in the phases of the Floquct eigenvalues (see, 
e.g., Fig. 1 of Ref. [26J|). For odd N, each partial CDT is 
associated with a perfect degeneracy of the phases of the 
Floquet eigenvalues while, for even N, some degeneracies 
(but not all) are incomplete provided Uq ^ 0. As in the 
case of modulated V , these incomplete degeneracies are 
caused by the existence of the Fock state \AN = 0). 
Consequently, partial CDT is generally more complete 
for odd iV than for even N. The Floquet spectrum in 
the case of large- amplitude modulation of U and weak 
interaction has been studied in depth in Refs. [22|, |26| . 
We refer to these references for further discussion. 

Finally, we point out that it is possible to create meso- 
scopic Schrodingcr's-cat-like states [NOON-like states 
[49J ] . i.e., states proportional to (\N) + e l9 \ — N)), where 
8 is a phase] using the state-dependent CDT. In this 
scheme, we assume that UqN/Jq -C 1 and choose \N) as 
the initial state. We modulate U at a frequency ui that 
corresponds to a partial CDT between |JV) and |iV — 2), 
that is, Jo [(Ui/u))(N - 1)] = 0. At this frequency the 
phases of the Floquet eigenstates ip^ ' , which are very 
close to NOON states, become degenerate [50]. By de- 
tuning from this partial CDT, we have a coherent oscilla- 
tion (with period T) between ijr N ' and ip% . As a result, 
the initial state \N) evolves into a NOON-like state at 
t = T(2n-l)/4, withn = 1,2,3,.... With increasing the 
absolute value of the detuning, the period T decreases 
but the amplitudes of the components other than | ± N) 
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FIG. 9: Time evolution of the normalized popula- 
tion imbalance (AN)/N and its variance uan/N = 
N^ 1 ^(AN 2 ) — (AN) 2 under large-amplitude modulation of 
U. Here TV = 21 [panels (a) and (b)] and TV = 51 [panels 
(c) and (d)], and the initial state is \AN = TV). In the case 
TV = 21 we have set Ui/Jo = 10 and lj/Jo — 83.85, and in the 
case TV — 51 we have set Ui/Jo = 4 and uj/Jo — 83.4. Other 
parameters are Uo — Ji — Vo = Vi = 0. A coherent oscilla- 
tion between |TV) and | — TV) is realized by slightly detuning 



(±) 



from a partial CDT between Floquet eigenstates tp 



increase, so that the oscillation between the NOON states 
is disturbed. Therefore, w (more precisely, U\/uS) should 
be optimized. In Fig. [9j we show the time evolution of 
the normalized population imbalance (AN)/N and its 
variance a AN /N = N- 1 V / (AN 2 )- (AN) 2 for TV = 21 
and TV = 51 as examples. Here (AN) = (tp\AN\ip) and 
(ATV 2 ) = (tp\(AN) 2 \tf,') with ATV = c\ci-c\c 2 . These are 
optimized cases with the amplitude of the wiggles in the 
oscillation of (AN)/N being < 0.05. When (ATV) = 0, 
can/N is almost equal to 1, which is the largest pos- 
sible value; this is a unique property of NOON states. 
Note that the oscillation periods arc comparable in the 
two cases: T/T = 211.3 and T/T = 367.4 for N = 21 
and N — 51, respectively. This shows that an advantage 
of the present scheme is that the optimized T does not 
increase exponentially with TV unlike the tunneling pe- 
riod of the higher-order co-tunneling in the self-trapping 
regime. This may be understood by the fact that the 
static part of the interaction strength Uq is very small 
(UqN/Jo <1). A disadvantage is that we need to know 
the number of particles exactly and to fine-tune U\/oj. 
This scheme can be used regardless of the value of TV if 
U N/J « 1. 



VI. CREATING A NOON STATE BY AN 
ADIABATIC SWEEP 

In this section we propose another scheme to create 
NOON-like states. This scheme uses an adiabatic sweep 
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FIG. 10: (Color online) Probabilities p a (t) = \(ipg\V(t))\ 2 
[blue (dark gray) lines], Ph{t) = |('(/'h|v['(i))| 2 [red (medium 
gray) lines], and p g + ph [green (light gray) lines] as a func- 
tion of time for two different values of the sweep rate a. Here 
TV = 5, Uo/Jo = 4, and A., = 0.5 (and V = V\ = Ui = 0). 
The dotted lines correspond to the analytical prediction ob- 
tained using Eq. (|30|) . 



of the modulation frequency. It enables us to obtain 
NOON-like states with TV < 10 particles starting from 
the ground state ip g of the timc-independent Hamilto- 
nian Hq. The basic idea is to create an avoided crossing 
between the Floquet eigenstate corresponding to ip g and 
the one corresponding to the NOON-like eigenstate tph 
by time-periodic modulation, which changes the geom- 
etry of the (quasi)energy space to be periodic. Here, 
we modulate the hopping parameter J and set the tilt 
V = 0. Since the phase 4>j of the modulation does 
not affect the result, we choose <j>j = for defmite- 
ness. The timc-independent part Hq of the Hamiltonian 
H(t) = H Q + H T ^ (t) is given by Eq. 1(15)1. while the time- 
dependent part H<r u (t) is 



Hr u (t) = -2J Aj sinuit S x 



(27) 




sufficiently large value, the first crossing occurs between 
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FIG. 11: (Color online) Asymptotic value p g of the transi- 
tion probability as a function of (a) the inverse sweep rate 
1/a and (b) the modulation amplitude Aj for TV = 5 and 
Uo/Jo = 4 (and Vo = Vi = Ui = 0). We set Aj = 0.5 in (a) 
and aT$ /tv 2 — 0.005 in (b). The circles show the numerical 
results, and the solid lines show the semi-analytic results ob- 
tained from the Landau-Zener formula (|30|) . The initial time 
ti of the time evolution is chosen such that u>(ti)/Jo = 29 in 
Eq. d28|. 



For even TV, the crossing used in the creation of the 
NOON state is the one between ^jr N ' and ipQ ■ For 
odd TV, it is the one between ip N and ip\ . We con- 
sider the regime UoN/Jq ^S> 1, where ip N is a NOON- 
like state. The ground state tp g of H corresponds to 
ip Q (even TV) or ip[ (odd TV), and the eigenvalue of 
Hq corresponding to ip g is denoted by E g . Similarly, 

the NOON-likc eigenstate iph of Hq corresponds to ip N . 
State iph has the highest energy among symmetric cigen- 
states of i/o, and its eigenencrgy is denoted by Eh- Be- 
cause H ~ UoN 2 ^> J N ~ Ht^ , the eigenstates of H 
are almost equal to the Floquet eigenstates except near 
the crossing points. Therefore, \(ip g \ipo )| 2 — 1 f° r even 

N, \(iPM +) )\ 2 - 1 for odd TV, and |(^|V4 +) >| 2 - 1- 
As discussed in Sec. IIV Al when u is decreased from a 



<.(+) 



odd TV 



and 



for even TV and between 4>\, T and <t>\ for 



■ iui uvcuiv auu ucuwccii i/y^ a,iiu y/j 

51| . Therefore, in principle, this scheme can be 
used without knowing precisely the total number of par- 



ticles. The avoided crossing between the phases 



d+) 



vo 

N 



or 



a(+) 






and 



is approximately at w r , 



E h 



E g . In the 



5 case discussed earlier, this crossing corresponds 
to the rightmost circle in Fig. EJb). 

Let us take ijj g as the initial state. If we sweep ui adi- 
abatically across the avoided crossing, ip g undergoes an 
almost perfect transition to iph- We consider a linear 
sweep of the form 



tj(t) = LU TC 



at, 



(28) 



where w ros is the location of the crossing and a is the 
sweep rate. The initial and final times of the sweep are 
denoted by t{ and tf, respectively. 

In the following calculations, we set TV = 5 and 
Uo/Jo = 4. The avoided crossing is at ui/Jo — 28. In 
Fig. [TU1 we show the time evolution of the probability 
p g (t) = \{ipg\&(t))\ 2 [blue (dark gray) lines] at which the 
system stays in the initial state ipg an d the probability 
Ph(t) = \(i/jh\^(t))\ 2 [red (medium gray) lines] at which 
the system undergoes a transition to the target state iph- 
Note that p g + ph shown by the green (light gray) lines 
in Fig. [10] is very close to unity throughout the calcu- 
lations (the deviation is within 0.1%), and the system 
is, to a very good approximation, restricted to the sub- 
space spanned by the two states. Thus the cross ing can 
be described by the Landau-Zener (LZ) model |52l - (55l ]. 
We denote the modulation period at the crossing point 
by T rcs = 27r/w rcs . The difference between the phases of 
the Floquet eigenvalues at u>(t) is A(f> = (cj) N ' — cpQ X J ) ~ 
— (Eh — T5 g )(T w — T rcs ). Here, we shift the phase dif- 
ference so that the crossing at uj ~ w rcs is passed at 
t = 0, in accordance with the standard expression of the 
LZ Hamiltonian. For the linear sweep of Eq. (f2"8")l . we 
get T u (t) = 2ir/uj(t) ~ (27r/wr e8 )(l + at/cu rcs ). Here we 
assume that at <C w ros . We obtain the quasienergy sep- 
aration AE corresponding to Acf> near the crossing as 



AE = 



Acj) 

± res 



at , 



(29) 



where we have approximated w r , 



E h - E„. The di- 



agonal matrix elements Hh and H g of the LZ Hamil- 
tonian are thus Hh, g = ±AE/2, where the upper sign 
corresponds to Hh and the lower one corresponds to 
H g . We found that the off-diagonal elements Hh g and 
Hgh = HI of the effective Hamiltonian are to a good 

approximation given by H hg = -JoAj(iph\S x \tp g )/V2. 
Consequently, the asymptotic value p g of the transition 



probability p g (t), p g = lim t 



,(t), is H 



Pg = CXP 



= cxp 



-2tt 



'fcfll 



|ft(if fc - H g , 



(30) 



In Fig. 1111 we show the probability p g as a function of 
the inverse sweep rate 1/a [Fig. fTTT a)] and the modula- 
tion amplitude Aj [Fig. flTT b)]. Since p g (t) and Ph(t) 
continue to oscillate around the asymptotic value un- 
til far after the crossing (see Fig. [TU]) . wc calculate p g 
by taking the time average of p g (£) after its oscillation 
amplitude becomes small and almost time independent. 
These results are shown by circles in Fig. 1111 Semiana- 
lytical results obtained from Eq. (|30[) are shown by the 
red solid lines. For the parameters used here (N = 5 
and Uq/Jq = 4), wc have E g /J a = 12.31, E h /J = 40.31, 

K^fc|&# fl >| = 9.697 x HT 2 , and w res « E h - E g = 
28.00Jo. The agreement between the semianalytical and 
numerical results is very good. 

Finally, we examine the experimental feasibility of this 
scheme. According to Eq. (|3U|) . to obtain a NOON-likc 
state, wc should satisfy the adiabaticity condition: 



nJgA^ h \S x \i; g )f 



>1. 



(31) 



In addition, the initial and the final frequencies should 
be outside the crossing region. Since the range of u char- 
acterizing the crossing region is comparable to the level 
separation A = 2\Hh g \ at the avoided crossing, the initial 
time ti and the final time t{ of the sweep have to satisfy 
\uj(t[j) — uj ics \ = a\t[j\ > 2\Hh g \. Also the Landau- 
Zener formula is valid under this condition. Taking into 
account the adiabaticity condition (|3T1) . this leads to the 
requirement 



|*i|,if > 



V2 



^ 4,1(^15*1^)1 



(32) 



As an example, let us estimate the time scale given by 
this equation by using the parameters used in the exper- 
iment of Ref. [53. In this experiment, the frequency of 
the pair tunneling is 4Jq/C/o — 550 Hz for Uq/Jo = 5; 
thus T ~ 0.72 ms. If Aj = 0.5, the right-hand side of 
Eq. d221) is 9 ms for N = 6, 40 ms for N = 7, and 214 ms 
for N = 8. Therefore, a NOON state with N <7 could 
be created within an experimentally accessible time, pro- 
vided the value of uj can be controlled with a sufficiently 
high accuracy. We note that, more generally, an upper 
limit for JV for this scheme to work is N ~ 10. Since the 
width of the peaks in the probability distribution (in the 
Fock space) of if) g and iph scales as ~ N 1 ^ 2 , a few times 
N 1 / 2 should be larger than N in order to have an over- 
lap between ip g and ?ph and to have a significant nonzero 
value of \(ip h \S x \tjj g )\. 



In the present scheme, the modulation of the hopping 
parameter works much better than the modulation of 
the tilt. This can be seen using perturbation theory. A 
straightforward calculation shows that for odd number 
of particles (iph\S x \ip g ) ~ {Jo/U )^ N ~ 3 ^ 2 and for even 
number of particles (ip h \S x \ip g ) ~ (J a /U ) {N ~ 2)/2 . In the 
same way, perturbation theory shows that (ip' h \S z \ip g ) ~ 

(Jo/C/o) (Ar - 1)/2 for odd N and (i>' h \S z \i> g ) ~ (J /U ) N / 2 
for even N. Here ip' h is the antisymmetric eigen- 
state of Ho with the highest energy. Wc sec that 

Wh\SM a )\ 2 /\{MSAi>g)\ 2 - (J0/U0) 2 , and consc- 
qucntly, the off-diagonal elements of the LZ Hamiltonian 
are much smaller when the tilt is modulated than when 
the tunneling is modulated. 



VII. CONCLUSIONS 

In this paper, we have considered a time-periodically 
modulated two-mode Bose-Hubbard model. We have dis- 
cussed three types of modulations, one where the tunnel- 
ing amplitude is modulated, another where the interac- 
tion strength is modulated, and a third where the energy 
difference between the modes (tilt) is modulated. Wc fo- 
cused mainly on the self-trapping regime, characterized 
by UqN 3> Jo j and assumed that the amplitude of the 
modulation is small. It is known that a modulation of 
the tunneling amplitude can lead to a drastic reduction of 
the tunneling period [28[ . We found that a similar effect 
can be induced by modulating the interaction strength 
or the energy difference between the modes. We have 
analyzed this phenomenon using Floquct theory as the 
main tool. We found that regardless of the modulated 
variable, the system has resonances at some modulation 
frequencies, corresponding to either greatly reduced or 
enhanced tunneling periods. To a good approximation, 
the locations of the resonances can be obtained with the 
help of the energy eigenvalues of the time-independent 
part of the Hamiltonian. Consequently, the locations of 
the resonances are almost independent of whether the 
tunneling, interaction, or tilt is modulated. 

We found numerically that if the tunneling amplitude 
or interaction strength is modulated, the system has a 
wide resonance; that is, the tunneling period is greatly 
reduced in a wide range of modulation frequencies. This 
resonance is present also in the case of modulated tilt, 
but it is much narrower. Furthermore, the behavior of 
the tunneling period as a function of the modulation fre- 
quency is not smooth in this case; see Fig.[5J These differ- 
ences can be explained using Floquct theory. The pres- 
ence of resonances is related to avoided crossings of the 
phases of the Floquet eigenvalues. In the case of a mod- 
ulated tunneling matrix element or interaction strength, 
the avoided crossings correspond to Floquet eigenstates 
with the same parity. In the case of modulated tilt, these 
avoided crossings correspond to eigenstates with opposite 
parity. In Sec. IIVC1 it is shown that due to this differ- 



ence, a wide smooth resonance cannot be obtained in the 
case of modulated tilt. 

We have also analyzed cases where the interaction en- 
ergy is weak in comparison with the tunneling energy, 
UqN/Jq < 1, and the modulation amplitude of cither 
the interaction strength or the tilt is large. Under these 
conditions, tunneling can be suppressed at some specific 
modulation frequencies. This phenomenon, the coherent 
destruction of tunneling, has been extensively studied in 
the literature. We concentrated on a property of CDT 
that has received less attention in the previous studies, 
namely, the Floquet spectrum of a system where the tilt 
is modulated. As expected, wc found that the suppres- 
sion of tunneling takes place when the phases of the Flo- 
quet eigenvalues become degenerate. For an even number 
of particles the suppression is more complete than that 
for an odd number of particles. 

Finally, we have proposed two ways to create a NOON 
state. One is based on coherent oscillation resulting from 
a detuning from a partial CDT caused by modulated in- 
teraction strength. An advantage of this method is that 
the tunneling period does not increase exponentially with 
the total number of particles N. The other method is 
based on sweeping the modulation frequency of the tun- 
neling term adiabatically. This scheme requires neither 
precise knowledge of the number of particles nor fine- 
tuning of the modulation frequency. We have shown that 
by using the latter method and the parameters of a recent 
experiment J56|, it is possible to obtain NOON states of 
N < 7 particles. 

It is known that the mean-field theory of the time- 
periodically modulated two-mode Bose-Hubbard model 
shows chaotic dynamics (e.g., Refs. 0, El E El H HH, 
[22,153). In the future, it would be interesting to study 
the connection between the Floquet spectrum of the orig- 
inal quantum system and the chaotic mean-field dynam- 
ics. Another interesting problem to study would be the 
quantum dynamics determined by a time-periodically 
modulated Hamiltonian in the presence of dissipation. In 
particular, the engineered dissipation leading to squeezed 
states proposed in Ref. [58[ is of interest. 
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Appendix A: Effective hopping parameter for 
modulated J 

Here we derive the effective tunneling amplitude in the 
limit of large-amplitude tilt modulation. The system fol- 
lows the Schrodingcr equation 



iiP(t) = H(t)i>(t), 



(Al) 



with 



H(t) = -2J S X + U S 2 Z + V{t)S z (A2) 



and V{t) given by Eq. ©. We go to a rotating system 
by defining 

$(t) = e ia M § *ip(t), (A3) 

where 

a(t) = f dr [Vo + Vi sm(ujT + fo)\ (A4) 

Jo 

Vi 
= Vot H [coscj>v — cos(ut + 4>v)\- (A5) 

Using this, the Schrodingcr equation becomes 

4{t) = H{t)i>{t), (A6) 

where 

H(t) = -2J (cos[a(t)} S x - sin[a(t)] S v ) + U S 2 Z . 

(A7) 

Assuming that the modulation period T u = 2ir /uj is the 
shortest time scale in the system, it is possible to obtain 
an effective Hamiltonian by averaging over T u as 



1 



H AVE (t) = — H(T)dT 

J-uj Jt 



t+T^ 



(A8) 



= -2Jf(t)S x - 2jf(t) S y + U Sl (A9) 



The effective tunneling amplitudes are defined as 

reSr.\ J 



J?(t) = 7T co S [a(r)}dT 

1 oo Jt 



Jf(t) 



Jo 



t+T u 



sin[a(r)] dr. 



(A10) 



(All) 



Instead of calculating J£ ff (£) and Jy B (t) separately, we 



write 



Jfit) - Uf(t) 



J e 



Vl 



IL 



t+T u 



^ T gi[VbT-^- cos(wt+0 v -)] _ 



(A12) 

This integral can be calculated easily using the equation 

00 
e * Z co S7 = J2 J n {z)e in ^+%\ (A13) 



where Jniz) arc Bessel functions of the first kind. We 
thus obtain 

jf{t)~ l Jf{t) 



r„ sin — 



v^ (V±\ e" l ( w * + ^'-5) V 



V 



u) J Vq + nui 



( Yl ]e l ^ co 



fee Z. 
(A14) 



In the special case Vb/w = k E Z, the original tunneling 
amplitudes J^ = Jo and J y — are replaced by effective 
ones, 



(^) 



Jf = Jo Jfe ( ^ ) 



where V\ is non-zero. 



fc (0V + 77) -cos<?V 



k [4>V + 2J cos^v 



(A15) 
(A16) 



Appendix B: Effective hopping term for modulated 

[/ 

In the case of large-amplitude modulation of the inter- 
action strength, the coherent destruction of tunneling is 
state-dependent [26[ . Here, we derive the effective Hamil- 
tonian for this case. 

We start from the time-dependent Schrodinger equa- 
tion (|A1[) with the Hamiltonian 



where U(t) is given by Eq. ([7]). For simplicity, we set 
V = 0. As in Appendix [Al we go to the rotating frame 
by defining 



^{t)=e la(t}§ h'(t), 



(B2) 



where 



a(t) = / dr[U + U x sin (lot + fa)] 
Jo 

= U(jt H [cos 4>u — cos (uit + </>u)] 

UJ 



(B3) 



Thus the Schrodinger equation becomes iip(t) 
H(t)ij){t), with 



H(t) = -Jo 



g ia(t)(2S,+l) + e -ia(t)(2S,+l)§_ 



, (B4) 



where S± = S x ± iS y . We have used the equations 
[SlS+] = S+{2S Z + 1), [SIS.] = -(2S Z + 1)£_, and 
S x = {S+ + S-)/2 to obtain 



1 



£«*(*) (2S.+1) + e -ia{t)(2Sz + l)g_ 



(B5) 



By time averaging over one modulation period T u , the 
effective Hamiltonian reads 






t+T^ 



H(t)cIt 



= - Jo[S+A + A^S. 



(B6) 



H(t) = -2 J S X + £/(t)S= 



(Bl) Here A is defined as 



A\AN) = < 



T u 



Jk 



■ sm 



^(Aiv + i; 



E ^ 



n— — 00 



— (AA + 1 



e i[£/o(t+5) + ^f cos </.[/] (AW+1) 



-f— - — IA7V), ^(AA + l^Z, 

U (AN + 1) -hw 1 " u 



U t 



(AN + 1) 



o i^L(A7V+l)cos0 [7 „-ife(0c/ + f ) 



|AA), 



f'„ 



(AA + 1) = k e Z, 



(B7) 



r 



where we have used the equation S Z \AN) = depends on AN. If (m — 2 1 -Have |to) = (here we assume 

(AA/2)| AA) and {\AN) ; AJV = -A, -A + 2, -A + m > without loss of generality), we get (m\H AVE \m - 

4, . . . , A} is the basis of the system. In this basis Have 2) = (-m + 2\H A ve\ - m) = (-m\H A v E \ - m + 2) = 0. 

is a tridiagonal matrix. Note that A, unlike Eq. (|A14[) . In the special case (Ua/uj)[(m — 2) + l] = k G Z, the con- 



dition for partial CDT between states \m) and \m — 2), 
(m — 2\HAVE\'m) = 0, can be written as 



Jk 



— (m — 1) 

UJ 



= 0. 



(B8) 



If this equation holds, the Hilbert space can be written 
as a direct sum of three uncoupled subspaces, spanned 

by{\N),\N-2),...,\m)}, {|m- 2), |m), . . . , | - m + 2)}, 
and {| - m), | - m - 2), . . . , | - N)}. 



[1] K. Hammerer, A. S. S0rensen, and E. S. Polzik, Rev. 

Mod. Phys. 82, 1041 (2010). 
[2] M. Saffman, T. G. Walker, and K. M0lmer, Rev. Mod. 

Phys. 82, 2313 (2010). 
[3] A. Cronin, J. Schmiedmayer, and D. E. Pritchard, Rev. 

Mod. Phys. 81, 1051 (2009). 
[4] K. Bongs and K. Sengstock, Rep. Prog. Phys. 67, 907 

(2004). 
[5] C. Lee, J. Huang, H. Deng, H. Dai, and J. Xu, Front. 

Phys. 7, 109 (2012). 
[6] D. H. Dunlap and V. M. Kenkre, Phys. Rev. B 34, 3625 

(1986). 
[7] F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, Phys. 

Rev. Lett. 67, 516 (1991). 
[8] M. Holthaus, Phys. Rev. Lett. 69, 351 (1992). 
[9] M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 
[10] F. Kh. Abdullaev and R. A. Kraenkel, Phys. Rev. A 62, 

023613 (2000). 
[11] M. Holthaus and S. Stenholm, Eur. Phys. J. B 20, 451 

(2001). 
[12] H. L. Haroutyunyan and G. Nienhuis, Phys. Rev. A 64, 

033424 (2001). 
[13] C. Lee, W. Hai, L. Shi, X. Zhu, and K. Gao, Phys. Rev. 

A 64, 053604 (2001). 
[14] G. L. Salmond, C. A. Holmes, and G. J. Milburn, Phys. 

Rev. A 65, 033623 (2002). 
[15] S. Kohler and F. Sols, New J. Phys. 5, 94 (2003). 
[16] H. L. Haroutyunyan and G. Nienhuis, Phys. Rev. A 70, 

063603 (2004). 
[17] A. Eckardt, C. Weiss, and M. Holthaus, Phys. Rev. Lett. 

95, 260404 (2005). 
[18] G.-F. Wang, L.-B. Fu, and J. Liu, Phys. Rev. A 73, 

013619 (2006). 
[19] C. E. Creffield, Phys. Rev. Lett. 99, 110501 (2007). 
[20] C. E. Creffield and F. Sols, Phys. Rev. Lett. 100, 250402 

(2008). 
[21] C. Weiss and N. Teichmann, Phys. Rev. Lett. 100, 

140408 (2008). 
[22] M. P. Strzys, E. M. Graefe, and H. J. Korsch, New J. 

Phys. 10, 013024 (2008). 
[23] X. Luo, Q. Xie, and B. Wu, Phys. Rev. A 77, 053601 

(2008). 
[24] D. Witthaut, F. Trimborn, and S. Wimberger, Phys. Rev. 

A 79, 033621 (2009). 
[25] J. Wang and J. Gong, Phys. Rev. Lett. 102, 244102 

(2009). 
[26] J. Gong, L. Morales-Molina, and P. Hanggi, Phys. Rev. 

Lett. 103, 133002 (2009). 
[27] Q. Xie and W. Hai, Phys. Rev. A 80, 053603 (2009). 
[28] G. Watanabe, Phys. Rev. A 81, 021604(R) (2010). 
[29] A. R. Kolovsky, Europhys. Lett. 93, 20003 (2011). 
[30] K. Kudo and T. S. Monteiro, Phys. Rev. A 83, 053627 

(2011). 
[31] G. Watanabe, S. Yoon, and F. Dalfovo, Phys. Rev. Lett. 

107, 270404 (2011). 



[32] I. Brouzos and P. Schmelcher, Phys. Rev. A 85, 033635 
(2012). 

[33] N. Gemelke, E. Sarajlic, Y. Bidel, S. Hong, and S. Chu, 
Phys. Rev. Lett. 95, 170404 (2005). 

[34] H. Lignier, C. Sias, D. Ciampini, Y. Singh, A. Zen- 
esini, O. Morsch, and E. Arimondo, Phys. Rev. Lett. 99, 
220403 (2007). 

[35] R. Gommers, V. Lebedev, M. Brown, and F. Renzoni, 
Phys. Rev. Lett. 100, 040603 (2008). 

[36] C. Sias, H. Lignier, Y. P. Singh, A. Zenesini, D. Ciampini, 
O. Morsch, and E. Arimondo, Phys. Rev. Lett. 100, 
040404 (2008). 

[37] A. Alberti, V. V. Ivanov, G. M. Tino, and G. Ferrari, 
Nat. Phys. 5, 547 (2009). 

[38] A. Eckardt, M. Holthaus, H. Lignier, A. Zenesini, D. 
Ciampini, O. Morsch, and E. Arimondo, Phys. Rev. A 
79, 013611 (2009). 

[39] A. Zenesini, H. Lignier, D. Ciampini, O. Morsch, and E. 
Arimondo, Phys. Rev. Lett. 102, 100403 (2009). 

[40] C. E. Creffield, F. Sols, D. Ciampini, O. Morsch, and E. 
Arimondo, Phys. Rev. A 82, 035601 (2010). 

[41] A. Alberti, G. Ferrari, V. V. Ivanov, M. L. Chiofalo, and 
G. M. Tino, New J. Phys. 12, 065037 (2010). 

[42] G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, 
Phys. Rev. A 55, 4318 (1997). 

[43] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, 
Phys. Rev. Lett. 79, 4950 (1997). 

[44] C. Chicone, Ordinary Differential Equations with Appli- 
cations (Springer, New York, 1999). 

[45] G. Watanabe and C. J. Pethick, Phys. Rev. A 76, 
021605(R) (2007). 

[46] As we shall see later, for odd N these N — 2 
crossing types consist of one type where </ > |ajv'I=jv-2 
crosses <p!& N \=N successively and N — 3 types where 
( t >( \AN\=N-4 >■■■ , <P\an\=i crosses (j>\^ Nl=N once. For 
even N, there is one crossing type with 4>\an\=n-2 
N - 4 types with each <p\tli\=N-4 , • ■ ■ , 0|a]v|=2 
one type with <^iv=o- In the case of modulated J (V), 
the phase of the Floquet eigenvalue <^Aiv=o undergoes 
avoided crossings with the phases of the eigenstates of 
even (odd) parity. 

[47] Based on the behavior of the phases of the Floquet eigen- 
values, there are resonances also at w/ Jo — 11 and 14. 
These are too narrow to be observed in Fig. [2] 

[48] For Uo — 0, the effective Hamiltonian [see Eq. (|A9f) ] 
vanishes when the condition for CDT is satisfied: 
Jo(Vi/w) — 0. When this equation holds, the phases of 
all the Floquet eigenvalues vanish, and the degeneracy is 
complete irrespective of whether N is even or odd. 

[49] H. Lee, P. Kok, and J. P. Dowling, J. Mod. Opt. 49, 2325 
(2002). 

[50] For even N this degeneracy involves three states if 



and 
and 



Uo — 0: Viv an( i another Floquet eigenstate whose 
quasienergy is identically zero. The third state disturbs 
the desired oscillation between ip N . The third state can 
be easily lifted to change the threefold degeneracy to a 
twofold one by introducing nonzero, but small, Uo such 
that UoN/Jo < 1. 
[51] This can be understood by the fact that, with decreasing 



U, (p 



(+) 



for even N and 



/,(+) 



for odd N deviate from zero 



n 

at the lowest rate and <j> N deviates at the highest rate 
among symmetric states. 
[52] L. D. Landau, Phys. Z. Sowjetunion 2, 46 (1932). 



[53] C. Zener, Proc. R. Soc. London, Ser. A 137, 696 (1932). 
[54] E. Majorana, Nuovo Cimento 9, 43 (1932). 
[55] E. C. G. Stiickelberg, Helv. Phys. Acta 5, 369 (1932). 
[56] S. Foiling, S. Trotzky, P. Cheinet, M. Feld, R. Saers, A. 

Widera, T. Miiller, and I. Bloch, Nature (London) 448, 

1029 (2007). 
[57] K. W. Mahmud, H. Perry, and W. P. Reinhardt, Phys. 

Rev. A 71, 023615 (2005). 
[58] G. Watanabe and H. Makela, Phys. Rev. A 85, 023604 

(2012). 



